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Abstract 

In two previous publicatio n] W » ^ , we have demonstrated that stationary rotation of magnetized 
plasma about a compact central object permits an enormous number of different MHD instabilities, 
with the well-known magneto-rotational instabilh y^liWii 5 ! as just one of them. We here concentrate 
on the new instabilities found that are driven by transonic transitions of the poloidal flow. A 
particularly promising class of instabilities, from the point of view of MHD turbulence in accretion 
disks, is the class of trans-slow Alfven continuum modes, that occur when the poloidal flow exceeds 
a critical value of the slow magnetosonic speed. When this happens, virtually every magnetic/flow 
surface of the disk becomes unstable with respect to highly localized modes of the continuous 
spectrum. The mode structures rotate, in turn, about the rotating disk. These structures lock 
and become explosively unstable when the mass of the central object is increased beyond a certain 
critical value. Their growth rates then become huge, of the order of the Alfven transit time. These 
instabilities appear to have all requisite properties to facilitate accretion flows across magnetic 
surfaces and jet formation. 
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I. INTRODUCTION 



In Fig. Q a Magnetized Accretion-Ejection Structure is shown which illustrates the prob- 
lem we wish to address in this paper, viz.: How does an accretion flow about a compact 
object first crosses the magnetic configuration and then turns the corner with respect to 
the accretion disk to produce jets? In ideal MHD, plasma and magnetic field stay together 
(frozen in field) , so that a sizeable resistivity is needed for the flow to detach from the mag- 
netic field. This involves anomalous dissipation. Hence, the basic problem is to find relevant 
local instabilities producing the necessary MHD turbulence. 




FIG. 1: (a) Stationary end state from simulation with YA( ' - (Versatile Advection Code): disk 
density surfaces (brown), jet magnetic surface (grey), helical field lines (yellow), accretion-ejection 
particle trajectory (red); (b) Accretion flow detaching from the disk. [From computations by Casse 
and Keppen ai 7 1ii 8 1 ]. 

Our model is shown in Fig. |21 An axisymmetric configuration of nested magnetic / flow 
surfaces with magnetic field indicated by the vectorial Alfven speed b and velocity v, having 
both toroidal and poloidal components, surrounds a compact object of mass M* in the origin. 
Note that the usual tokamak configuration is obtained for M* = 0, whereas accretion disk 
geometries may have flat (thin disk) as well as round (thick disk) poloidal cross-sections.^ 
This model considers laboratory and astrophysical toroidal plasmas on an equal footing by 
exploiting the scale independence^ of the MHD equations. 

In order to obtain a stationary equilibrium situation, we assume that the accretion flow 
speed is much smaller than both rotation speeds of the disk. We then need to determine the 
stationary equilibrium flows (Sec. 2) and, next, the local instabilities driven by the transonic 
flow (Sec. 3). We analyze this problem from two angles: 

(a) Asymptotic analysis for small inverse aspect ratio (e <C 1); 

(b) Large-scale exact numerical computations of the equilibria and instabilities. 

These will be discussed in reverse order since the numerical results suggest the relevant 
approximations that may be made. 
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FIG. 2: Transonically rotating magnetized disk about compact object. [From Ref. I[2]l ]. 

A difficulty encountered is that transonic transitions upset the standard equilibrium- 
stability split. We will discuss this in Sec. 2 under the heading of transonic enigma. 
The gravitational parameter exploited here is defined as follows: 



for parallel flow 



Rv% 



R M 2 B 2 

This parameter measures the deviation from Keplerian flow (for which T = 1) 



(1) 



II. TRANSONIC EQUILIBRIUM FLOWS 

We exploit a variational principle for the stationary axisymmetric equilibria 

determining the poloidal flux tp and the poloidal Alfven Mach number squared M 2 ( = 
pVp/B 2 y This involves five arbitrary scaled flux functions Aj(^>) : 

x'{ip) = pvp/ B p , (derivative poloidal stream function) (2) 



H(if>) = hvZr= ; + — L ---$R l n i - , =, (Bernoulli function) (3) 

y 1 2 p B 2 7 - 1 p 2 y^WTZ 2 1 J K 1 

S(ip) = p~ 7 p, (entropy) (4) 

= R^ 1 [v v — {x'/p) By] , (—derivative electric potential) (5) 

K(ip) = R \v v — (l/x')B v \ , (pol. vorticity & current density stream function) (6) 

which have to be fixed by whatever observational evidence is available. Note that the 
flux function (JHJ) has been renamed K (instead of the usual L) to eliminate the confusion 
derived from the frequent occurrence of the misnomer 'specific angular momentum' in the 
astrophysics literature (see references in Ref. [ 2.1). This is important since one of the essential 



problems in accretion disk dynamics is precisely the transport of angular momentum. 
The stationary states are then determined by minimization of the following Lagrangian: 

sfcdV^O, £ = ^(l-M 2 )|Wf -^7--^+ 113 2 , (7) 
J 2R 2K n Yl M -?M 2 ^ 1-M 2 ' v; 

where Uj(Ai(ijj)] R, Z) are simple algebraic expressions. The Euler equations provide the 
solutions i/j(R, Z) and M 2 (R, Z) of the core varia bles. [A generalization of this variational 



principle to two-fluid plasmas is given in Ref. 
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The transonic enigma mentioned above is due to the fact that the flows suddenly 
change character from elliptic to hyperbolic at the transonic transitions. As a result, standard 
(tokamak) equilibrium solvers diverge in the hyperbolic regimes! We circumvent this problem 
by calculating in elliptic regimes beyond the first hyperbolic one. Obviously, the payoff is 
that we cannot approach the transonic transitions but have to infer what has happened 
there from the changes in the dynamics found in the 'transonic' elliptic regimes. 

The pleasing side of the transonic enigma is that the time- dependence of the linear waves 
and the spatial dependence of the nonlinear stationary states are intimately related. This is 
seen by comparing the wave spectra, which cluster at the slow, Alfven, and fast continuum 
frequencies {cuf}, {^}, oJp = oo for highly localized modes (Fig. EJ), and the corresponding 
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FIG. 3: Cluster spectra of the waves. 



slow, Alfven, and fast hyperbolic flow regimes delimited by critical values of the square of 
the poloidal Alfven Mach number (Fig. EJ): 
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FIG. 4: Flow regimes of the stationary states. 



By means of the transonic equilibrium solver FINESSE (described in Ref. \12]) typical 
equilibria in the first trans-slow elliptic regime have been computed for tokamak and 
accretion disk (Figs. and |HJ) for a representative choice of the flux function parameters: 




FIG. 5: Equilibrium for tokamak (r = 0). 



Note that the accretion disk equilibrium, in contrast to the tokamak, has the density 
peaking on the outside to produce overall equilibrium on the flux / flow surfaces with respect 
to the gravitational pull of the compact object in the center. For the sake of the spectral 
calculations (Sec. 3), the two equilibria have been chosen such that the safety factor q{ip) is 
a monotonically increasing function. 
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III. LOCAL TRANSONIC INSTABILITIES 



With equilibrium flows, the overall spectral structure of MHD waves and instabilities 
is determined by the split in forward and backward waves so that the local waves cluster at 
the Doppler shifted continuous spectra: 



±cus + k • v , 



±u; a + k ■ v 



±00 



(8) 



These are embedded in a monotonic spectral structure for ID equilibria, as shown in Fig. [7| 
A long-standin g puz zle about the nature of the singular frequencies fin (= k ■ v) has been 
clarified in Ref. I[13]t In the Eulerian description, these frequencies give rise to the Eulerian 
entropy continua (£Ie), not perturbing the pressure, velocity, and magnetic field, and, hence, 
absent in the Lagrangian description. 
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FIG. 7: Schematic spectral structure for ID stationary equilibria. 

To solve for the local transonic instabilities, we exploit the Frieman-Rotenberg equation, 
F stat ic(£) + V ■ [p(y ■ Vv)£] + p(u + zv ■ V) 2 | = . (9) 

The linear operator is no longer self-adjoint, so that overstable modes occur, in particular in 
2D axisymmetric equilibria through coupling of the poloidal modes e m ^. 

In order to solve for the transonic continuum modes we exploit localization on sepa- 
rate magnetic / flow surfaces, 



This gives rise to an eigenvalue problem for each surface, 



(10) 



A V = B V 



(11) 
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where the matrices are defined by 



( R 2 B 2 . _ 2 .£?2 
-F-g/-F-(M 2 -M c 2 )- 



fi 2 



z(M 2 - M 2 C ) 



B 2 
P 2 







pRB, 
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Tp 
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%P T{M 2 -M^ 
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8 fpRB z 

\ B 2 



TM 2 C B 2 T + p 



/ R 2 B 2 

1 Jpu-TM)—gf(jpu-MT) -ta^pu 



B EE 



d({M 2 - M 2 c )yd P 



\ 



;i2) 



V 



ia^fpui 



(13) 



(^pu - FM)B 2 (^pu - MT) ) 



and the Doppler shifted frequency uj = uj — nQ in a frame rotating with Q ( ^ v v /R !) . 

The overall result of the analysis of this eigenvalue problem (111)1 is that the continuum 
modes are always unstable in the trans-slow ( M 2 > M 2 ) flow regime! This is due to the 
poloidal derivatives indicated by the terms [d (. . .)] multiplied with factors (M 2 — M 2 ). The 
instability mechanism involves coupling of the Alfven and slow continuum modes. 

Fig. |H1 shows the full complex spectrum of trans-slow Alfven continuum 'eigenvalues' for 
a Tokamak (r = 0) equilibrium with the radial coordinate s = ip 1 ^ 2 as a parameter. The 
colors indicate on which magnetic / flow surface the modes are localized. The modes rotate 
clockwise for Reu > and anti-clockwise for Re a) < 0, and have growth rates of the order 
of a few percent of the inverse Alfven transit time. 
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FIG. 8: Complex eigenvalues of n = 1 transonic instabilities for tokamak (r = 0). 

The counterpart for an accretion disk (r = 2) equilibrium is shown in Fig. El Now 
the rotating continuum modes lock to give purely exponential growing modes (Re a; = 0) 
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over a sizeable range of magnetic / flow surfaces. Their growth rates are huge, in the order 
of ten to twenty percent of the inverse Alfven transit time! Consequently, these modes have 
enough time to saturate during a finite number of revolutions of the plasma. 
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FIG. 9: Complex eigenvalues of transonic n = 1 instabilities for accretion disk (T = 2). 



Analysis of the dispersion equation by small inverse aspect ratio expansion (e C 1) 
is suggested by the numerical results, which exhibit dominant coupling of the six Alfven 
and slow continuum modes S^-i) ^mi ^m+i around the e = degeneracies at the rational 
surfaces q = —m/n. This analysis confirms that the trans-slow Alfvenic continuum modes 
are unstable at, or close to, the rational surfaces for all toroidal mode numbers n. 

For a very massive central object (r ^> 1), the growth rate in the limit |n| , \m\ — > oo 
becomes 



to 



2p 



MT 



(14) 



which is far in excess of the Alfven frequency. 

Since these modes are localized both radially (because of the continuous spectrum) and in 
the angle $ (and ip) they are perfectly suitable to produce the turbulence that is needed along 
the accretion flow at the inner edge of the accretion disk with respect to the central object to 
detach the flow from the magnetic field. 

IV. CONCLUSIONS 



(1) In the presence of poloidal rotation, the singular structure of the MHD continua transfers 
to the equilibrium, so that linear waves and non-linear stationary equilibrium flows are no 
longer separate issues. 
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(2) Complete spectra of waves and instabilities have been computed for tokamaks and ac- 
cretion disks with transonic flows exploiting our new computational tools FINESSE (for 
equilibria) and PHOENIX (for stability). 

(3) We have found a large class of instabilities of the continuous spectra of transonic ax- 
isymmetric equilibria for M 2 > M 2 . 

(4) These instabilities may cause strong MHD turbulence and associated anomalous dissi- 
pation, breaking the co-moving constraint of plasma and magnetic field and facilitating both 
accretion and ejection of jets from accretion disks. 
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